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It is argued why thermodynamic approximations in terms of dressed propagators may, at larger coupling 
strength, be better behaved than perturbative results, and why in hot QCD the hard thermal loop 
approximation of the thermodynamic potential cannot be expected to work close to the phase transition. 

1. Introduction 

The thermodynamic potential is a quantity of central importance in statistical field theory, and a lot 
of interest has been devoted to its calculation. At high temperature, the perturbative expansion of the 
thermodynamic potential is known to 5th order in the coupling for the <j> A theory jjj , , for QED || , and 
for QCD ||. However, for physically relevant questions the applicability of these results is limited; 
they are meaningful only when the coupling is so weak that the thermodynamic potential is very close 
to its free limit. Already for moderate values of the coupling the accuracy of the approximation does not 
improve by the inclusion of higher order terms. 

This feature of the perturbative expansion is similar to the behavior of asymptotic series. Consider, 
e.g., the integral Z(g 2 ) = dx exp{— L}, where L = | x 2 + g 2 x 4 : which may be seen as a caricature 
of the partition function Z in the </> 4 theory ||. Z(g 2 ) can be expressed in terms of a modified Bessel 
function. In the spirit of perturbation theory, however, one would expand the quartic 'interaction part' 
of the exponential, and perform the integrals before the sum to obtain an expansion in powers of g 2 . 
Being related to the number of diagrams in the perturbative calculation of the partition function in the 
</> field theory, the coefficients of this series, (-l) n /n! fdxx in exp{-x 2 /2} ~ (-l) n n!, alternate in sign 
and grow asymptotically like the factorial of the order n. Therefore, by the ratio test, the radius of 
convergence is zero. Nonetheless, the coefficients contain relevant information. In fact, the accuracy of 
the approximation of Z(g 2 ) by the truncated series does improve by higher-order terms as long as the 
truncation order is smaller than an optimal n* ~ 1/g 2 . 

The divergence of the perturbative expansion of Z(g 2 ) is obvious already in its definition: in the 
complex g 2 plane, Z has a cut along the negative axis. Due to this essential feature one may presume 
that some of the properties of Z(g 2 ) and its expansion translate, at least qualitatively, to path integrals 
and, hence, to the thermodynamic potential £1 = — T\nZ in field theories. For example, the analysis 
of the (non-) convergence of the perturbative expansion of f2 demonstrates how the optimal truncation 
order decreases with increasing coupling. It is noted that for QCD close to the transition temperature 
already the next-to leading order perturbative contribution cannot be reconciled with the lattice data. 

To motivate the following calculation of the thermodynamic potential in terms of the propagator, I 
want to point out another feature of the above toy model. Firstly, one can realize that the analog of the 
propagator, d(g 2 ) — (x 2 )/Z, where the weight of the average is exp{— L}, has also a cut along the negative 
g 2 axis. The same holds for the the 'self-energy' p(g 2 ) = d^ 1 (g 2 ) — c? _1 (0). This suggests to consider Z 
as a complex function of p. It turns out that the mapping p — > g 2 — > Z unfolds the discontinuity along 
the negative g 2 axis. This analytic aspect is a formal motivation to choose, in the physics context, a 
representation of the thermodynamic potential in terms of the self-energy, or the dressed propagator, as a 
starting point to derive approximations. This is, of course, also supported by physical intuition. Dressed 
Greens functions are expected to be a preferable basis for nonperturbative approximations, since they 
resum entire subsets of 'bare' graphs. In particular, as known in other fields of physics, the excitations 
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described by the full propagator may be considered as quasiparticles whose effective properties encode 
already a part of the interaction effects. The idea to understand the thermodynamics of a strongly coupled 
system in terms of its quasiparticles was successfully demonstrated for hot QCD in phenomenological 
models where the self-energies are approximated by their asymptotic values, see, e. g., |?]]. 

In the following, a more systematic approach towards a quasiparticle description of the thermodynamics 
of field theories is given, starting from ^-derivable approximations. To outline this concept, the scalar 
theory is considered first in the next Section, before then turning to gauge theories. 



2. ^-derivable approximations, scalar theory 

The thermodynamic potential was first expressed as a functional of the exact propagator for a non- 
relativistic fermion system [|J. This representation is the starting point for the so-called conserving 
approximation scheme J9|, in which diagrams are resummed in thermodynamically consistent subsets. In 
the scalar theory, the thermodynamic potential can be expressed by the functional 
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evaluated at its stationary point, — 0. This condition is equivalent to II/2 = 5&/5A. It defines 

the full propagator A, which is related to the self-energy II by Dyson's equation, A -1 = Aq 1 — II. $ is 
given by the sum of the dressed 2-particle irreducible bubble diagrams; for the interaction <7q/4! 4> 4 , 
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from which the expansion of II/2 is obtained diagrammatically by cutting one line in each graph. The 
loop truncation of these skeleton expansions yields consistent ('^-derivable') approximations 

At leading-loop order, the self-energy is just a mass term^. The ^-derivable approximation can be 
obtained explicitly and it is equivalent to the large- N limit of the 0(iV)-symmetric scalar theory 
fl2| . Considering, for simplicity, the massless theory, and working in 3 — 2e dimensions in the MS scheme, 
the self-energy is to be determined selfconsistently by 



n« = 12Q=|f-^[--ln?| + ll + / where uj — (k 2 + II/;) 1 / 2 
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The right hand side contains a divergent contribution ~ n«/e. This term, which implicitly depends on the 
temperature T, is the vacuum part of the resummed quantum fluctuations, while the medium contribution 
is finite due to the thermal distribution function rib. The divergence is absorbed in the bare coupling 
constant. The renormalized coupling g is defined as the scattering amplitude at a momentum scale s. 

Resumming the chains^, X = X + 12 XIX ' yields # 2 ( s ) = 9o - l9W( s )[ l / e ~' [n {- s /P 2 ) + 2 ]/(16tt 2 ). 
Expressing now go by g renders eq. (0) finite and independent of fi, 
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where the coupling g has been defined at — s = T 2 . This gap equation has a nontrivial structure. For g 
smaller than g w 10, there are two solutions. The smaller one is related to the perturbative result which, 
by expanding the right hand side in n«/T 2 , is reproduced up to order g 3 . The second solution, related 
to the tachyon is for small g 2 exponentially large and thus of no physical relevance. As g — > g, the 
two solutions approach the same value, II = Tlu(g) ss 4T 2 , and the leading-loop approximation can no 
longer be justified. For even larger values of the coupling, (Q) has no solution. 

2 Therefore, the resulting approximation is similar to that in screened perturbation theory juj with the 'appropriate' choice 
for the mass parameter, introduced there to reorganize the perturbative expansion. 

s This yields only 1/3 to the first coefficient of the /3 function, which is an artefact of the truncation of the loop expansion. 
The two crossed scattering diagrams have to be omitted here since they would induce graphs beyond the superdaisy subset. 
This is consistent with the fact that they are suppressed in the 0(N — > oo) theory, where the leading-loop approximation 
becomes exact. 
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At leading-loop order, the ^-contribution <E>m = i ifilluAu combines with the first term in eq. (0), 



n b (co/T) n 



2 



, : ,J 1287T 2 ■ ^ 



I emphasize the relative factor of 4 between the two terms under the sum-integral; it ensures that the 
divergent T-dependent contributions of these terms cancel: (11^ /87r) 2 [— l/e+| 2/e] = 0. The first integral 

in eq. (||) is the thermodynamic potential of a system of quasiparticles with mass II*/ 2 . Expanded in 
II;z/T 2 ~ g 2 , it over- includes the perturbative 0(g 2 ) term by a factor of two. This is compensated by the 
second integral - in fact, the approximation (|^) reproduces the 0(g 3 ) result. The term ~ IL^ is a remnant 
of the resummed vacuum fluctuations. It leads to a minimum in the pressure, p = — ft/V, considered as a 
function of the self-energy, see Fig. 1. Due to the stationary property of ft«, the location of this minimum 
coincides with II;;, and hence indicates also where the leading-loop approximation breaks down. 

This observation is interesting. With regard to the next section, assume we had solved the gap equation 
only perturbatively, to leading order: 11/; ss II* = (g T) 2 /2A (without resummation, the coupling remains 
bare). This is equivalent to the hard thermal loop (HTL) approximation of II. Approximating now the 
^-contribution by $* = j $IT A* yields the same functional form (§) for ft* (IT) as in the selfconsistent 
approximation^. Therefore, the maximal value of the self-energy, and hence for the coupling, beyond 
which the approximation cannot be meaningful (which could not be surmised from II* = (goT) 2 /24) can 
be inferred also from the HTL approximation ft* (II*), which is important for a later argument. 

Another useful thermodynamic quantity is the entropy density, s — dp/dT. In the leading-loop $- 
derivable approximation, it can alternatively be obtained from a 1-loop functional of the dressed propa- 
gator JljJ, E3, which is not specific to the </> 4 theory. In this case, the result is equivalent to the entropy 
of a system of quasiparticles with mass II; 1 / 2 in the selfconsistent approximation, or (II*) 1 / 2 in the HTL 
approximation. Therefore, the entropy is a monotonously decreasing function of the quasiparticle mass, 
see Fig. 1. It shows, in contrast to the pressure, no sign of the breakdown of the approximation. This 
observation is consistent with the thermodynamic relation, since to reconstruct p from s one needs the 
T-dependence of the self-energy. It is also plausible on general grounds. While the pressure and the 
self-energy contain vacuum contributions, which become important at larger coupling, the entropy is 
determined only by thermal excitations. 



3. Gauge theories, QCD 

Aside from the technical complications to solve a set of coupled nonlocal Dyson equations selfconsis- 
tently, it may not be desired in gauge theories. In the <I>-derivable approximation scheme, 'bare' diagrams 
of all orders in the coupling are grouped in subsets to preserve thermodynamic consistency. Since this 
procedure distinguishes the 2-point functions in the hierarchy of Greens functions, it cannot be expected 
to also preserve gauge invariance (unless solved exactly to all orders in the loop expansion). Another, 
technical, issue is how to accomplish the renormalization of resummed self-energies at finite tempera- 
ture. Approximating the self-energies at leading order in the coupling makes the renormalization trivial. 
This would correspond to the considerations at the end of the Sec. 2, and lead to similar consequences. 
For example, only the leading order of the pressure could be reproduced upon expanding the resulting 
expression, while the (gauge dependent) higher order terms would not match the perturbative ones. 

However, it turns out that approximating the self-energies by their HTL contributions, which are gauge 
invariant, suffices to resum the leading contributions. When calculating the thermodynamic potential, 
this additional approximation requires to account for the ^-contribution^ more carefully than in the cj> A 
theory. For simplicity of the argument, consider first the case of QED, 

ft = iTr[ln(-D- 1 ) +L>n] - Tr[ln(-,S , - ;L ) + 5S] -*[£>, 5], n = 2 S<S> / '5D , E = -S$/6S, (6) 

4 Of course, it is not the same function of the coupling, and it can reproduce the perturbative result only up to leading order 
(as an aside, the next-to leading order term comes out too small by a factor of ^). Still, it is a useful approximation since 

the quasiparticle mass is a measurable quantity. Note also that an ansatz <E>* = 3<jig(^ A*) 2 has to be abandoned since it 
leads to uncompensated T-dependent divergenc es in the thermodynamic potential. 

5 This is not an issue in the entropy calculation |1 4| because the leading loop entropy-functional is independent of <&. 
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where D = (D 1 — II) 1 is the photon propagator, and S = (S 1 — E) 1 is the electron propagator. The 
Lorentz and spinor indices are implicit here, and summed over in the traces. At leading loop order, 




Both in the selfconsistent and in the perturbative approximation, the ^-contribution to the thermody- 
namic potential could be evaluated by closing the external legsQ of, either, the photon or the electron 
self-energy, or a linear combination of those traces. This is not the case in the HTL approximation, where 
terms of order of the external momentum squared are neglected for the contributions ~ T 2 . However, 
when traced over all momenta, these terms contribute to the same order, e 2 T 4 , as the HTL self-energies. 
This does not indicate that the HTL approximation is inutile hereQ, but rather suggests the solution of 
this riddle, namely by keeping track of all relevant terms in At order e 2 , is given by a double 
sum- integral over an expression with a numerator N = K 2 — Q\ — Q 2 , where K is the photon momen- 
tum, and Qi,2 are the electron momenta. One of the three terms in N would be neglected in the HTL 
approximation, e.g., closing the legs of the photon HTL self-energy H* amounts to drop K 2 . Summing 
over the cases where one of the terms in N is dropped yields 2N . Hence, $* = 4 Trn*£)* — \ TrE^S 1 * is 
equivalent to <!>;; up to terms of higher order. The resulting approximation p"5| for f2, 

fi* = iTr[ln(-L>*- 1 ) + in*L» i ] - Tr[ ln(-5* ~ x ) + \ E*5* ] , (8) 

is in close analogy (notice the factors i in the trace terms) to the expression in the scalar theory. 

For QCD, a similar (though due to the non-Abelian topology slightly more involved) analysis of the 
HTL ^-contribution was given in |15| |. Let me here argue directly for a SU(iV c ) gluon plasma (including 
ghosts to compensate the unphysical degrees of freedom) why to expect the expression 

= \ Tr c [ ln(-£>* _1 ) + \ n*Z>* ] - Tr c [ ln(-Dj ostfi )] . (9) 

The transverse and longitudinal components of the gluon HTL self-energy Tl*(po,p) are given by di- 
mensionless functions of po/p, times the asymptotic gluon mass squared, M 2 = \N c g 2 T 2 . Therefore, 
M 2 d M 2 ln(— D*^ 1 ) = — H*D*. Hence, the factor i between the two terms in the gluon trace is indispens- 
able to cancel their individually divergent contributions, which are ~ Mg/e n on dimensional grounds. 




Figure 1. Pressure and entropy in the HTL approximation, as functions of the thermal mass scale in the cf> 4 
theory (left), and for SU(iV c ) (right). Shown also are the leading and next-to leading perturbative results. 
In the 4 theory, the selfconsistent and the HTL approximation have the same functional dependence on 
the respective self-energy. In the hatched region the approximations are not justified. 

6 This is the same strategy as in the (p 4 theory, see footnote^. Evaluating, instead, the functional <E>;j with the perturbative 
propagators would also in the present case result in uncompensated T-dependent divergences. 

7 It is noted that the leading contribution in the final result for the thermodynamic potential comes from large momenta 
close to the light cone, where the HTL approximation is justified. 
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The result (H) has a striking similarity to the scalar case. Its expansion reproduces the 0{g 2 ) result, but 
under-estimates the 0(g 3 ) term by a factor of j. The pressure p* has a minimum at Mg sa 3T 2 (see 
Fig. 1). Similar to the (j) 4 theory, the dip is caused by a term ~ M 4 implicit in ([)]). It originates from the 
vacuum contributions, which for larger coupling become as important as the thermal contributions. This 
is in accordance with the general argumentation in Sec. 2, and not specific to the HTL approximation. 
I therefore conjecture also for QCD, where the solution structure of the selfconsistent Dyson equation 
cannot be studied explicitly, that the leading- loop approximation may not be meaningful for M g M g . 
Again, this could not be inferred from the HTL entropy TL JT^| which decreases monotonously with 
Mg. This conjecture is supported by lattice data. At M g , p* is approximately 80% of the free pressure, 



Figure 2. The pressure as a function of the temper- 
ature in the SU(3) theory. The shaded band shows 
the lattice data p6[ , Jl7| . The hatched band for p* 
represents the variation of the result when determin- 
ing the temperature dependence of Mg by, either, 
the 1- or the 2-loop running coupling in the MS 
scheme, assuming T c = l.lAAjjg, and varying the 
renormalization scale 1 < £i/(itT) < 4 as in p4| . In 
accordance with the conjecture in the text (but not 
shown here), also the HTL entropy TL (T) starts to 
deviate clearly from the lattice data below 3T C Jl4| . 

a value reached by the lattice data at 3T C (T c is the confinement temperature). At this temperature, 
indeed, p* g {T) starts to match the data, see Fig. 2, as does the entropy s g 1TL |Q. 

In conclusion, it was shown that the HTL approximations of thermodynamic observables has a con- 
siderably improved range of applicability in the large coupling regime compared to perturbative results. 
For hot QCD, the approximation works down to temperatures of a few times above T c . This encourages 
developments to extend, similar to the approach |lS|], the formalism to finite chemical potential, where 
no lattice data are available so far. 
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